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Abstract -Up to now there has been no rehable method to calculate the Casimir force when 
surface roughness becomes comparable with the separation between bodies. Statistical analysis 
of rough Au films demonstrates rare peaks with heights considerably larger than the root-mean- 
square (rms) roughness. These peaks define the minimal distance between rough surfaces and 
can be described with extreme value statistics. We show that the contributions of high peaks to 
the force can be calculated independently of each other while the contribution of normal rough- 
ness can be evaluated perturbatively beyond the proximity force approximation. The developed 
method allows a reliable force estimation for short separations. Our model explains the strong 
hitherto unexplained deviation from the normal Casimir scaling observed experimentally at short 
separations. 
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' Introduction. — The Casimir force [T] attracts in- 
qreasing attention novifadays since modern technology al- 
lows dimension control at distances < 100 nm where this 
force becomes operative (see [11[3] for a review). Indeed, 
modern micro/nano-clcctromechanical (MEM/NEM) en- 
gineering is now being conducted at the micron to nanome- 
ter scale and has attracted interest in the Casimir force 
[3]. MEM devices such as vibration sensors and switches 
are now routinely made with parts a few micrometers in 
^ize, and have the right size for the Casimir force to play 
a role. This is because MEM systems have surface areas 
lp,rge enough but gaps small enough, for the force to draw 
components together and possibly lock them permanently 
-\ an effect known as stiction. Such permanent adhesion (in 
addition to capillary adhesion due to the water layer) is a 
common cause of malfunctioning of MEM devices [SHT] . 

In this range of separations the force appears mainly 
due to quantum fluctuations of the electromagnetic field 
(zero-point field) in the interacting bodies while at larger 
distances classical (thermal) fluctuations become increas- 
ingly important [SKH]- The famous Casimir formula Fc = 
{tt^ / 240) {he/ d^) gives the force (per unit area) at tem- 
perature T = between two ideally reflecting semi-spaces 
separated by the distance d. The force measured in recent 



experiments (see [3] for a review) can deviate significantly 
from the ideal case because the temperature is finite, the 
bodies are not ideal reflectors, and the distance between 
them is not well defined. Considerable efforts were made 
to improve the Casimir formula. Indeed, the more de- 
tailed description is based on the Lifshitz formula, which 
accounts for actual optical properties of interacting bodies 
and nonzero temperature. The optical data were included 
in the calculational procedure [IHIIII] • Although the ther- 
mal correction to the Casimir force is rather controversial 
[121 US], it is not important for the short distances dis- 
cussed in this paper. 

An important correction to the Casimir force that is 
not accounted for by the Lifshitz formula is the roughness 
correction. The surfaces of real bodies are rough, which 
makes the distance between them not well defined. The 
first attempts to account for roughness [2] were based on 
the proximity force approximation (PEA) . In this approx- 
imation the real surfaces are replaced by flat patches and 
the force was calculated as the sum of forces between op- 
posite patches, treating such pairs as parallel plates. Eor 
the dispersive forces the PEA was applied for the first time 
by Derjaguin |15lll6j . The approximation is justified when 
the separation d is much smaller than the local curvature 
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radius and size of patches. It is well suited for smooth large 
bodies, but works worse for roughness corrections. It was 
noted [T7| that in order to apply the PFA to rough bodies 
the roughness correlation length ^ (typical features size on 
the surface) must be larger than the separation, ^ 3> c?. 
Then the result found in |14| will be true for small root 
mean square (rms) roughness, w <^ d. However, in most of 
the experimental situations the condition ^ 3> d is broken 
and more elaborate theory has to be used to calculate the 
roughness correction. This theory was developed in refs. 
[18,19 . It treats the roughness contribution through sec- 
ond order perturbation theory in w/d. The theory showed 
larger corrections than those predicted within the PFA. In 
fact, the correction is very important at short separations 
and has to be carefully included for interpretation of the 
force experiments exploring short distance ranges. 

The Casimir forces between a gold covered sphere and 
plates of different roughness were measured for separa- 
tions from 20 to 200 nm [50]. The films with larger 
rms roughness at short separations demonstrate signif- 
icant (more than 100%) deviations from the theoretical 
expectations based on the perturbative roughness correc- 
tion. Empirically it was established that the minimal dis- 
tance between two rough bodies (distance upon contact) 
is c?o ~ i.7{w + Wsph), where Wgph is the sphere's rms 
roughness. Because do is the minimum separation dis- 
tance, the perturbative correction must be smaller than 
K{w + WsphY /d'^ = 0.07i^, where K ^ 10 is a large nu- 
merical factor (due to sharp behavior of the force with 
the distance). It was concluded [5D] that at short sepa- 
rations the perturbation theory fails. These experimental 
results still did not get a theoretical explanation. More- 
over, we are facing a problem: there is no reliable method 
to estimate the roughness correction when d becomes com- 
parable with the rms roughness w. In this paper we pro- 
pose a method to address this problem by combining the 
PFA and perturbation theory approaches. Although we 
will prove the applicability of this method specifically for 
gold films, we believe that similar approaches can be de- 
veloped for other materials, after detailed analysis of the 
roughness statistics obtained, e.g., in terms of scanning 
probe microscopy techniques. 

Statistics of rough surfaces. — The distance upon 
contact do was discussed in detail for gold films [H]. The 
films deposited with different thicknesses have different 
rms roughnesses due to kinetic roughening processes. For 
all these films atomic force microscope (AFM) images were 
recorded for large area (of up to 40 x 40 ^m^) with lateral 
resolution of 4-10 nm. This information allows a detailed 
analysis of the roughness statistics. The probability to 
find a height of a local feature smaller than some value z 
can be presented in a general form 



P{z) = 1 



(^) 



(1) 



describes the roughness distribution in a convenient way, 
which makes it possible to calculate the contributions from 
peaks and troughs as will be shown later. It was already 
noted [21] that the cumulative distribution P{z) for gold 
films cannot be described satisfactorily by any known dis- 
tribution at all z but asymptotically at large z it can be 
fitted with generalized extreme value distributions |22] . 

We performed a special analysis of the AFM surface 
data presented in ref. [21] to reveal the best asymptotic 
distribution at large \z\. In this limit the phase 4>{z) is 
much more convenient for analysis than P{z). This is 
because P{z) approaches very fast or 1 in the limit \z\ — > 
oo. Indeed, we can present the phase as 



r) = - In [1 - P{z)] 



(2) 



where P{z) is extracted directly from the images. The 
function (f>{z) for an 1600 nm thick gold film is shown in 
fig. [TJ The inset shows the probability density function 
f{z) ~ dP/dz = {1— P)d(j) / dz . Similar behavior is realized 
for all investigated gold films. It is clear that for large 
positive z the logarithm of the phase can be fitted with a 
linear function 



ln0(z) = A + Bz, 



(3) 



and similarly for large negative z. With this (f> the prob- 
ability to find a feature larger than z behaves asymptoti- 
cally as a double exponential 



1 — P{z) ^ exp 



' exp 



(4) 



where for convenience we introduced the "phase" 0(z) as 
nonnegative and nondecreasing function of z. The phase 



where 13 and ^ are the scale and location parameters re- 
spectively. This behavior is a characteristic feature of the 
Gumbel distribution [33] , which is an example of extreme 
value statistics. In this paper only gold films were an- 
alyzed and therefore we cannot draw conclusions on the 
roughness statistics of other materials. However, the ex- 
treme character of the statistics allows us to hope that 
this behavior is more general. 

Roughness correction to the Casimir force. 

We can imagine a rough surface as a large number of as- 
perities with heights w and lateral size ^, and occasional 
high peaks and deep troughs. These peaks (or troughs) are 
high in the sense that their height is considerably larger 
than w, say > 3w. The situation can be visualized as a 
lawn covered with grass and occasional high trees standing 
here and there. In this paper we propose a method to cal- 
culate the roughness correction to the Casimir force on the 
basis of this separation. Namely, the asperities with the 
height ^ w can be taken into account using perturbation 
theory, without the use of the proximity force approxima- 
tion. On the other hand, for high peaks the local distance 
between interacting bodies becomes considerably smaller 
and one cannot use perturbation theory anymore. Because 
high peaks are rare the average distance I between them is 
large. If this distance is so large that I ^ d, then we can 
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Fig. 1: (Color online) The "phase" as a function of z for a 1600 
nm gold film. The open circles are the actual data extracted 
from the AFM image using eq. ((2]). At large positive and large 
negative heights logj^g (?i(2) is well fitted by linear functions of 
z as is shown by the straight lines. The curved line is a poly- 
nomial fit at intermediate z. The inset shows the probability 
density function /(z). It demonstrates significant deviation 
from a normal distribution. 

calculate the contribution of these peaks independently 
of each other, as it is assumed in the PFA. However, this 
contribution has to be calculated beyond the perturbation 
theory. It has to be stressed that the interaction of a sep- 
arate peak with a flat surface can be taken into account 
precisely using developed numerical or analytical methods 

i- 

The number of asperities N with the height d\ > Sw 
and lateral size ^ on the area is given by the equation 

m 



N 



r2 



The average distance between these peaks is 



L 

N 



(di)/2 



(5) 



(6) 



In order to fulfill the condition of PFA applicability / 3> d, 
we can choose the parameter di from the interval < 
di < do, where c?o is the maximal peak on the area L^. 
The best choice for di will be discussed later. Similarly, 
one can introduce the average distance I' between deep 
troughs 



I' 



(7) 



where d'l has to be chosen from the interval 3w < d[ < d'^ 
to fulfill the condition V ^ d and d[f is the deepest trough 
on the area L^. 

Here we consider the general case where we are in- 
terested in the Casimir force between two plates with 
rough surfaces. As was explained in ref. this is 
equivalent to the interaction of a smooth plate with a 



rough one, which has the combined roughness topogra- 
phy h{x,y) = hi{x,y) + h2{x,y), where hi^2{x,y) are the 
topographies of the interacting plates 1 or 2. Therefore all 
the equations above have to be applied to the combined 
roughness profile h{x,y). 

Let us assume for a moment that the PFA can be applied 
to any roughness topography. Then the force between the 
plates can be calculated using the standard definition of 
the averaged function 



da -<ii di 

/..,+ /.,.+ / 

di —d'n —d'. 



dzf{z)F{d- z), (8) 



where f{z) is the probability density function, and wc sep- 
arated high peaks (first integral) , deep troughs (second in- 
tegral), and the normal roughness contribution (third inte- 
gral) . For the moment we do not specify the force between 
the interacting patches F{d) separated by the distance d. 
The last term can be calculated using the perturbation 
expansion F{d - z) = F{d) - F'{d)z + F"{d)z^/2\ + ... 
and we find for this term 



Fid) / dzfiz) + 



F"{d) 
2! 



dzf{z)z\ (9) 



-d\ 



The first and second integral on the right are 1 and uP' , 
respectively, if one extends the integration limits to infin- 
ity. When the applicability of the PFA breaks down the 
second term in ^ (with infinite limits) can be generalized 
as follows [TTI 



FpT{d) 



F"{d) 
2! 



d^k 
{2-kY 



p{kd)a{k). 



(10) 



Here <y{k) is the Fourier spectrum of the roughness correla- 
tion function. The function p{kd) measures the deviation 
from the PFA. When the PFA is applicable this function 
is p{kd) ~ 1 and we reproduce eq. Outside of the 

PFA applicability we can use for p{kd) expressions found 
in [HllTg. 

The term J-pxid) is the roughness contribution to the 
force treated as a perturbation theory correction. As we 
know already the contribution of high peaks (or deep 
troughs) may not be accounted for by the perturbation 
theory. However, in this case we can account for the 
peaks (troughs) independently and the contribution can 
be presented as the first (second) term in ([5]). Taking 
into account the change of the integration limits we find 
the contributions due to high peaks, TppAid), and due to 
deep troughs, J-^'ppj^{d), 



TpFA{d) 



do 



Jdzfiz) F{d-z)-F{d) + F'{d)z 



F"{d) .,2 
2! 



(11) 
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T'pp^{d) 



J dz,f{z) F{d- z)- F{d)+F'{d)z- ^^z^ .(12) 



The final expression for the force that includes the total 
roughness contribution can be presented as 

T{d) ^ Fid) + FpT{d) + FpFAid) + F'ppAid). (13) 

Here F{d) is the force between flat surfaces and the other 
three terms are the different roughness corrections. 

The same force F{d) is used to calculate FpFA{d) and 
Fppj^{d), which implies that high peaks are described as 
pillars with flat faces. However, this approximation is 
not necessary. If the peaks can be considered as inde- 
pendent, then the interaction of each peak with the fiat 
surface can be described precisely (numerically) or approx- 
imately with an appropriate force F{d) in eqs. pip and 
([T2|) , taking into account the actual geometry of the peak. 
For example, high peaks can be considered as pillars with 
spherical caps of radius ^/2. As we will see below for the 
description of the experiment |20] it is sufficient to use the 
simplest model for the peaks (flat faces). 

At this point an important question is: with what pre- 
cision can we calculate the roughness corrections? The 
term —F'"{d)z^/3\, which is neglected in the Taylor ex- 
pansion of F(d — z), allows an estimation of the error in 
FpT- In the distance range that we are interested in here, 
20 < d < 100 nm, the force F{d) behaves with the dis- 
tance as F{d) = A/d"- , where A is a constant and a « 3.5 
1241. Then the error is estimated as 



AFpT = 7 



a{a + l){a + 2) 
3! 



/"w 



F{d), (14) 



where 7 is the skewness of the distribution f{z). The 
data shown in fig. [T]give the largest 7 ~ 1.285 among the 
investigated films and we estimate the error as AJ^pt ~ 
18. 55 {w/d)^F{d). The latter means that the perturbation 
theory correction has meaning at least for d/A > w. The 
minimal distance between rough surfaces do depends on 
the area of nominal contact L^, but even for L as small 
as 1 /im^ the condition do/w > A is usually fulfilled [?T] . 
Therefore, we can now draw the important conclusion that 
the perturbation theory correction (|10p has meaning up to 
the point of contact between interacting rough surfaces. 

The precision with which we calculate the contribution 
of the high peaks is defined by the condition of applicabil- 
ity of the PFA to these peaks. This condition is l{di) ^ d 
and we have for the error 



AFpFA = {d/l)FpFA- 



(15) 



As we already mentioned we have to choose di such that 
the condition /(di) ^ d is true and, therefore, the correc- 
tion (|lip makes sense. Similarly, we can define the error 
AF'ppji^ for the contribution of deep troughs (fT^ . The rel- 
ative error in (jlip increases with the distance, but we have 



to keep in mind that Fp fa decreases very fast with d and 
the absolute error stays small. The parameters di and d'l 
can be chosen rather arbitrarily if the conditions l{di) 3> d 
and l{d'i) » d are fulfilled. A practical recipe could be 
di = max [iw, {do + 'w)/2] and d'l = max [3w, (dp + w)/2]. 
It has to be noted that the contribution of deep troughs 
is always small, but we keep it for the sake of generality. 

Results. — The roughness corrections (jlOp - (fT^ were 
deduced for the force between two rough parallel plates. In 
most of the experimental configurations the sphere-plate 
geometry is used. We can find the result for this configura- 
tion if the sphere's radius is large, i? S> d. This condition 
is typically true when the roughness effect is appreciable 
and we can apply the PFA to the total force F{d). The 
same equations (|10p - (|13p can be applied but now we have 
to understand F{d) as the force between a smooth sphere 
and a smooth plate, approximated by 



F{d)^2'nRE{d) R:^d, 



(16) 



where E{d) is the Casimir-Lifshitz energy per unit area 
for the parallel plate configuration [5]. We neglect the 
thermal effect (T = 0) at short separations [12] . However, 
we use measured optical properties of gold films [TT] to 
account for the actual material properties. 

We evaluated the force and all the roughness correc- 
tions to compare it with the experimental data [20]. The 
Lifshitz force F{d) was calculated for the sphere radius 
i? = 50 /xm using the optical data for sample 3 in |llj . 
The roughness effect was estimated for 800, 1200, and 1600 
nm Au films and a Au covered sphere. Here we present 
the results for the 1600 nm film. The roughness char- 
acteristics for combined sphere-plate AFM images were 
presented in ref. [H]. For the 1600 nm film they are: rms 
roughness w = 10.1 nm, the correlation length ^ = 42 nm, 
and the distance upon contact do = 50.8 ± 1.3 nm. The 
last value was determined by electrostatic calibration [20] . 
It is preferable to use this value, because dg determined 
from the roughness topography has a larger uncertainty 
[H]. We used di = {w + do)/2 = 30.5 nm. According 
to eq. ^ it corresponds to the average distance between 
high peaks, I ~ 380 nm. Note that the effective area of in- 
teraction is L^, with L = 2100 nm [31]. For deep troughs 
the calculation details are less important. For the given 
L we found dg = 24.6 nm. Since dg < 3w the troughs are 
not deep enough and can be taken into account perturba- 
tively. Therefore, in this specific case there is no need to 
introduce J'pfa- 

The results are presented in fig. [21 One can see that the 
solid (blue) line, which shows the result of our approach, is 
in agreement with the experimental data within the exper- 
imental errors. This is in contrast with the perturbation 
theory approach that failed to explain the data [20] . This 
is demonstrated in the inset, which shows different com- 
ponents of this force. At short distances the contribution 
of high peaks (2, red) is so large that it dominates the 
whole force. In this case a few peaks become very close 



p-4 



Roughness correction to the Casimir force beyond perturbation theory 






50 60 70 80 90 100 

d (nm) 

Fig. 2: (Color online) The force between a Au covered sphere 
{R = 50 ^m) and a plate (1600 nm thick Au). The open 
(green) circles are the experimental data from [5D]. The vertical 
and horizontal bars show the experimental errors for a few 
points. The solid (blue) line is the result of our model. Naive 
application of the PFA to the force between rough bodies based 
on eqs. ^ and (|16|) is shown by the dashed (red) curve. The 
inset shows different components of the force. 1 (black) is the 
force F{d) between smooth surfaces calculated according to the 
Lifshitz formula. 2 (red) is the contribution of the high peaks 
according to eq. 3 (blue) is the perturbation theory 

correction according to eq. (|10p . The sum of all three curves 
gives the solid line in the main panel. 

to the opposite body so that the force diverges at d — )- do. 
There can be very few high peaks but their contribution 
cannot be neglected. On the other hand the contribution 
of high peaks disappears very fast when the distance be- 
comes larger. 

We used two different models to calculate the contri- 
bution of high peaks in eq. (|lip . In the first model the 
peak was considered as a pillar with a flat face. In the 
second model the peak had a spherical cap of radius ^/2. 
The interaction of the cap with a plate was taken into 
account according to ref. [IS], where the proximity force 
approximation is not used. We found a negligible differ- 
ence between the two models of peaks. The reason is the 
following: When the distance d — rfo ^ then the contri- 
bution of the peaks is very small due to their small area of 
interaction. When d approaches do or d — do <?C the con- 
tribution of high peaks becomes significant, but the shape 
of the peaks is not important anymore, because the PFA 
is valid in this limit. 

Naive application of the proximity force approximation 
according to eqs. (O and gives the dashed line (red) 
in fig. [2j It is interesting to note that this line is also 
in agreement with the experimental data. At the shortest 
separations both curves coincide, because the dominating 
high peaks can be treated with the PFA. At larger dis- 
tances the perturbative contribution becomes important 
and the PFA result lies below the solid line as it should 



be [T7]. However, the difference between these two curves 
is within the experimental errors. Perturbation theory ac- 
counts for the non-additivity of the Casimir force, whereas 
the PFA assumes it is additive. So this difference pro- 
vides an indication of the effect of the non-additivity in 
the roughness correction. It can be concluded that within 
the experimental error the experiment in ref. [20| was not 
sensitive to this non-additivity. 

Conclusions. — In conclusion, we developed a reliable 
method to include the effect of roughness of interacting 
bodies in the Casimir force at short distances when per- 
turbation theory fails. It was established that roughness of 
gold films can be described asymptotically (for high peaks 
or deep troughs) by extreme value statistics. In this case 
the rough surface can be presented as a large number of 
asperities with heights of the order of the rms roughness 
and a few occasional peaks, which are much higher than 
the rms roughness. The distance between high peaks is 
large so that one can calculate their contribution for each 
peak separately (using the PFA). The smaller asperities 
can be calculated using perturbation theory beyond the 
PFA. The contribution of high peaks is extremely impor- 
tant for short separations, where it dominates not only the 
perturbative roughness correction but also the force as a 
whole. Therefore, our result is interesting not only for 
the Casimir force but also for the problem of adhesion be- 
tween surfaces in general [26j , including wet environments 

We repeat that the method presented here solves the sig- 
nificant discrepancy between measurements of the Casimir 
force at short separations, and the results of perturbation 
theory [ID]. 
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